Charge and critical density of strange quark matter 
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The electric charge of strange quark matter is of vital importance to experiments. A recent in- 
vestigation shows that strangelets are most likely highly negatively charged, rather than slightly 
positively charged as previously believed. Our present study indicates that negative charges can 
indeed lower the critical density, and thus be favorable to the experimental searches in heavy ion 
collisions. However, too much negative charges can make it impossible to maintain flavor equilib- 
rium. {S0556-2813(99)02506-6} 

PACS numbers: 24.85.+P, 12.39.-x, 12.38.Mh, 25.75.-q 
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The last 20 years have witnessed an ever-increasing in- 
terest in strange quark matter (SQM) jD. Early in 1970's 
PL it had already been known that the appearance of 
s fraction in two flavor quark droplets provides an addi- 
tional degree of freedom, and thus could lower the energy 
of the system. In 1984, Witten conjectured that quark 
matter with strangeness per baryon of order unity might 
be bound || , which has aroused a great focus of investi- 
gations into the stability and detectability of SQM. In the 
framework of MIT bag model, Jaffe et al. find that SQM 
is absolutely stable around the normal nuclear density for 
a wide range of parameters |Q while other authors find a 
few destability factors . On application of the quark 
mass-density-dependent model, Chakrabarty et al. car- 
ried out a lot of investigations with significantly different 
results . However, the investigations by Benvenuto 
and Lugones give similar results to those in the bag model 
H . A more recent investigation indicates a link of SQM 
to the study of quark condensates pi. 

The most important way of the territorial searches for 
SQM is in the altra-relativistic heavy ion collisions [ pd| . 
Of vital importance to the experimental searches is the 
electric charge of SQM. Previously, it is generally be- 
lieved that SQM is slightly positively charged [Q. Con- 
trary to previous findings, a recent investigation shows 
that metastable strangelets are most likely highly nega- 
tively charged [Q . This significant result may, if really 
the case, have a serious impact on the high sensitivity 
searches in heavy ion experiments at the AGS and CERN 
facilities. Our present study, in a very simple manner, in- 
dicates that proper negative charges can indeed lower the 
critical density under which SQM can not maintain its 
flavor equilibrium, and thus no longer exist. 

Following Ref . Q , We assume the SQM to be a Fermi 
gas mixture of it, d, s quarks and electrons with chemical 
equilibrium maintained by the weak interactions: d, s <-> 
u + e + V e , s + u <-> u + d. 

At zero temperature, the thermodynamic potential 
density of the species i is 
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where i — u,d,s,e (ignore the contribution from neu- 
trinos); gi is the degeneracy factor with values 6 and 2, 
respectively for quarks and electrons. 

The corresponding particle number density is 
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For a given baryon number density n b and electric 
charge density Q, the chemical potentials fj, u , fi d , Msj an d 
ix e are determined by the following equations 



fJ-d = Ms = M> 
Mu + Me = Ml 

\{n u +n d + n s ) = n b , 

~3^ s ^ e 

The last two equations are equivalent to 

n u - n e — n b + Q, 
n d + n s + n e = 2n b - Q. 

We therefore define a function of fi e 

F(fj, e ) = n 2 (n d + n s + n e ) - Tr 2 (2n b - Q) 
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where 



M = Me + \Jm 2 u + [n 2 (n b + Q) + 1(^ 2 - m 2 ) 3 / 2 ] 2 / 3 . 
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Because m s > m u ^, the equation F(/j, e ) = for /i e has 
solution if and only if 

fi > m s , (12) 
F(p e ) < 0. (13) 

At the critical density n c of baryon number, the equal- 
ity signs in the above two inequalities should be taken. 
So we can easily find the equation 
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which determines how the critical density depends on 
charge density, namely, n c — f(Q). 

At n c , the strangeness fraction becomes zero. When 
the density decreases further, the equation group (^-||) 
which determines the configuration of the system has no 
solution. This indicates that n c is the lowest density for 
the possible existence of SQM. 

It is obviously difficult to get an explicit expression for 
the function f(Q). However, the inverse function f~ l (n c ) 
can be easily obtained from Eq. (pTj): 
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where 
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In order to make the inverse square root meaningful, 
we must require 



2\3/2 
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When taking the equality sign, and then solving the 
corresponding equation, we obtain the minimum critical 
density n cm i n . Substituting into Eq. (Bjh gives 
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where m! s and m' u are, respectively, the masses of s quarks 
and u quarks corresponding to the density n cm i n . 

If the charge density Q is more negative than Q m i n , 
there is no critical density corresponding to it. Therefore, 
Q m i n is the most negative charge density of SQM. 

To perform a concrete calculation of the critical density 
corresponding to Q > Q m in, we must stop here for a little 
while to discuss the consideration of quark confinement. 

Let's schematically write the QCD Hamiltonian den- 
sity as 



Hqcd = H k + qM q + Hi, 



(19) 



where H k is the kinetic term, Hi is the strong interaction 
terms between quarks, Mq is the current mass matrix 
with diagonal elements m u o, m^o, and m S Q. 

In the mass-density-dependent model, the strong in- 
teraction between quarks is mimicked by the proper de- 
pendence of quark masses on density. It is equivalent to 
expressing the Hamiltonian density as 



-ffequiv = H k + qMq 
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which is the same with that of free quarks except that 
the current mass matrix Mo is replaced by an equivalent 
mass matrix M with diagonal elements m q (q — u, d, s). 
From the equality 



(jib l-f^cquiv \rib) = (n&|-ffQCD|nb) , 
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where \nj,) is the state vector with baryon number density 
rib, we have 



(n b \Hi\n b ) 



J2( n b\Q(l\ n b) 

q 



m q0 + mi. 



(22) 
(23) 



The physical meaning of Eq. (g3|) is clear: the first 
term on the right is the original mass, i.e., the quark 
current mass, while the second term m\ is the inter- 
acting part mimicking the strong interaction between 
quarks. Obviously, m\ is flavor-independent and density- 
dependent. Because the characteristic of strong inter- 
actions is the confinement of quarks, we should require 
X\m nh ^mi - co. A popularly used parametrization 
form is @-|] 
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where B is a fixed constant determined by stability ar- 
guments. 

Because we are interested in the true ground state 
of SQM which might pass into a detector, we must re- 
quire, at zero pressure, the energy per baryon is less than 
Nh6 Fe /56 = 930 MeV for SQM and greater than 930 
MeV for two-flavor quark matter (where Maep e is the 
mass of 56 Fe) M. A little later, we will see that the 
negative charge density is not permitted to shift too far 
away from zero. We thus use the same method with that 
in Ref. H to estimate the absolute stability range for 
the paremeter B (the authors of Ref. M rename it C), 
which gives (69, 110) MeV-fm~ 3 . Taking the modest val- 
ues B = 80 MeV-fm -3 , m u0 = m d0 = 7.5 MeV, and 
to s o = 100 MeV, we can calculate the critical density by 
substituting Eq. (fj) into Eq. Q or @. The result 
is plotted in Fig. 1 with a solid line. It is obvious that 
n c is an increasing function of Q. Negative charges cor- 
respond to smaller critical densities. However, when Q 
decreases further toward negative direction, there will be 
no n c corresponding to it because SQM can not main- 
tain its flavor equilibrium at that case. Therefore, we 
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conclude that proper negative charges can lower the crit- 
ical density, and thus be favorable to the experimental 
searches. 

To explore if the above conclusion depends on mass 
formulas used, we now adopt the following expression: 

D . . 

m q = m q0 + — , (25) 

obtained for reasons of symmetry (particle-antiparticle 
conjugation, Lorentz transformations). Eq. ( p5| ) is obvi- 
ously in consistency with the properties of the mass (in- 
variant under particle-antiparticle conjugation, Lorentz- 
scalar; here the density is replaced by square of the four- 
current The constant D should be in the range 
of (113-130 MeV) 7 according to the method in Ref. ||. 
Taking D = (120 MeV) 7 and the same values of m 9 o with 
the above, we plot the corresponding critical density as 
a function of the charge density in Fig. 1 with a dashed 
line. 

We also perform a calculation in the bag model where 
the quark confinement is mimicked by adding to the en- 
ergy density expression an extra constant which is inter- 
preted as the vacuum pressure, while the quark masses 
are density-independent: 

m u , d = 7.5 MeV, (26) 
m s = 100 MeV. (27) 

The corresponding result is given in Fig. 1 with a dotted 
line. 

The shapes of the three lines in Fig. 1 are obviously 
similar to each other and so lead to the same conclusion 
as mentioned in the above. The significant difference 
is that the line corresponding to bag model calculation 
is much lower, or in other words, the charge density of 
stable bulk SQM in the bag model is nearly not allowed 
to be negative. This is due to the complete ignorance of 
the strong interactions in the bag model calculation. At 
high densities, they asymptotically become identical as 
expected. 

It should be pointed out that our method of reason- 
ing is different from what given by Schaffner-Bielich et 
al. Jl^ | . There the investigation was concerned with how 
small metastable strangelets (including explicit finite size 
corrections) look like and might decay for different scales 
of lifetime. None of the candidates found was being abso- 
lutely stable, but this fact, of course, depends on the bag 
parameter being employed. There it was not assumed, as 
done here, that the strangelets are in weak chemical equi- 
librium, but it was investigated how the metastable can- 
didates might look like if they are assumed to be stable 
against strong hadronic decay and subsequently against 
weak hadronic decay. These decay modes drive the more 
stable candidates to negative charges. Here we have as- 
sumed that the pieces of SQM are in perfect weak chemi- 
cal equilibrium and considered as having the possibility of 
absolute stability. Within the standard MIT approach, 



these pieces of stable bulk SQM would have a slightly 
positive charge. It is shown here, with the quark mass- 
density-dependent model, that the overall charge in fact 
might be slightly negative. Naturally, the present treat- 
ment is very simple, so the concrete values should not be 
taken very seriously, and further studies are needed. 

In conclusion, the present investigation is in favor of 
the viewpoint that proper negative charges can lower the 
critical density and consequently be favorable to the ex- 
perimental searches for strange quark matter in ultra- 
relativistic heavy ion collisions. 
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FIG. 1. Critical density n c vs electric charge density Q. 
The three lines correspond to different mass formulae (for 
details, see text). We can see that n c is an increasing function 
of Q, i.e., negative charges can lower the critical density. The 
points marked with "O" correspond to the most negative 
charge density permitted. 
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